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EQUILIBRIUM  AMO  STABILITY  PROPERTIES  OF  THE  SOLENOIOAL  LENS  BETATRON 


THoms  P.  Hughes  and  Brendan  B.  Godfrey 
Mission  Research  Corporation,  1720  Randolph  Road,  S.E. 
Albuquerque,  New  Mexico  87106 


Summary 

The  solenoldal  lens  betatnwF uses  a  series  of  solen¬ 
oids!  lenses  arranged  around  a  race-track  shaped 
drift-tube  to  provide  strong  transverse  focusing  for 
a  high-current  electron  ring.  Equilibrium  behavior 
of  the  circulating  beam  Is  examined. .for  parameters 
close  to  those  of  the  University  of  Mm  Mexico  mach¬ 
ine  currently  under  construction.  %The  tolerance  of 
the  beam  to  mismatches  In  the  toroidal  and  vertical 
fields  Is  evaluated  analytically  and  using  a  particle 
simulation  code.  The  linear  and  nonlinear  develop¬ 
ment  of  the  negative-mass  Instability  In  the  device 
Is  also  studied.  -.Stability  behavior  comparable  to 
that  In  a  conventional  betatron  (l.e.,  one  with  no 
toroidal  magnetic  field)  Is  found.  Growth  rates  are 
compared  to  those  obtained  from  an  analytic  model. 

Equilibrium  Properties 


In  order  to  confine  and  accelerate  high  currents 


In  a  betatron,  the  conventional  weak  focusing  must  be 
supplemented.  In  the  solenoldal  lens  betatron  (SIB)1 


this  Is  accomplished  through  the  use  of  periodic  sol¬ 
enoldal  lenses,  as  shown  In  Fig.  1.  The  SLB  thus 
differs  from  the  *mod1f1ed"  betatron,2  for  example, 
which  uses  a  uniform  toroidal  field.  For  a  matched 
equilibrium,  the  beam  In  the  SLB  reverses  Its 
pololdal  rotation  at  each  lens,  so  that  the  net 
pololdal  rotation  Is  zero.  This  requires  that  the 
beam  be  Injected  from  a  cathode  which  Is  shielded 
from  magnetic  flux,  so  that  the  beam  produced  has  no 
canonical  pololdal  angular  momentum.  In  the  region 
between  lenses,  the  beam  rotates  at  the  Larmor  fre¬ 
quency,  l.e.,  Qe/2y,  where  fi«  Is  the  nonrelatlv- 
Istlc  cyclotron  frequency  In  the  solenoldal  field, 
and  y  Is  the  relativistic  factor.  If  the  beam  emir, 
tance  Is  negligible,  then  force-balance  requires 


1  -  ",  ♦  P2'«  •  0 


where  ns  -  2vR2/y,|2aa,  Is  the  self-field  Index,  v 
Is  Budker's  parameter  (beam  current  divided  by  17 
kA),  R  and  a  are  major  and  minor  beam  radii,  respect¬ 
ively,  e  Is  the  beam  velocity  normalized  to  c,  p  Is 
the  ratio  Bg/Bz,  where  Bz  Is  the  vertical 
betatron  field,  and  tfe  Is  the  solenoldal  (toroidal) 
field  between  the  lenses.  We  have  assumed  that  the 
vertical  field  Index  Is  1/2.  From  Eq.  (1),  we  can 
compute  the  solenoldal  field  needed  for  a  matched 
beam,  given  the  other  quantities.  Sample  parameters 
for  Initial  experiments  at  the  University  of  Hew 
Mexico  JUNM)  are:  I)  •  1  m,  a  •  2  cm,  y  »  3,  v  ■ 
5.9*10"*  (100  Amps).  From  these,  we  compute  that 
Bg  •  80  Gauss  Is  required.  This  result  Is  expected 
to  be  accurate  In  the  limit  where  the  thickness  of 
the  magnetic  cusps  Is  much  less  than  the  distance 
between  them.  A  scenario  closer  to  that  envisioned 
for  the  um  device  Is  to  have  21  cm  long  solenoids 
spaced  10.5  cm  apart.  There  are  thus  20  solenoids 
evenly  spaced  around  the  1  m  major  radius  torus.  We 
model  the  magnetic  fields  by  using  the  exact,  numer¬ 
ically  computed  values  for  finite  length  straight 
solenoids.  One  "cell"  of  the  periodic  field  Is  shown 


■W5FI  Supported  by  the  Office  of  Naval  Research. 


In  Fig.  2.  Simulations  of  beam  behavior  In  this 
device  are  carried  out  using  the  electromagnetic  par¬ 
ticle  code  IVORY.*  We  find  that  to  minimize  the 
envelope  fluctuations,  the  field  In  the  middle  of 
each  solenoid  must  be  about  100  Gauss,  the  discrep¬ 
ancy  between  this  and  the  analytic  result  Is  presum¬ 
ably  due  to  the  thickness  of  the  cusps,  «Ai1ch  lowers 
the  average  field  value. 
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Figure  1.  Section  of  Solenoldal  Lens  Betatron.  The 
pololdal  velocity  Vpol  reverses  from  one 
solenoid  to  the  next. 
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Figure  2. 


Vector  plot  of  one-half  period  of  toroidal 
field,  showing  location  of  solenoldal 
windings. 


Mismatch  In  Solenoldal  Field 


To  check  the  tolerance  of  the  beam  to  mismatch 
In  the  solenoldal  field,  we  run  two  cases  In  which 
the  value  of  Bg  Is  doubled  and  halved  respectively. 

In  each  case,  we  ensure  that  the  beam  has  zero  canon-  oil 
leal  pololdal  angular  momentum.  The  resulting  envel- Jp 
ope  oscillation  amplitudes  are  shown  In  Fig.  3.  Them 

beam  oscillates  between  radii  of  2  cm  and  0.5  cm  for _ 

Bg  ■  200  Gauss,  and  between  2  cm  and  3  cm  for  Bg  , 

■  50  Gauss.  Experimentally,  It  should  not  prove  dif¬ 
ficult  to  avoid  these  large  mismatches.  Lty 


if 


* 


HI 


1:1 


lty  Codeajjj 
[Avail  and/or  »* 

I  Special  £ 


.*  190. 


94. 

>9.49  >9.99  9.99  9.99  9.49 

»(•») 

Figure  3.  fcxplltude  of  envelope  oscillations  for  (1) 
Be  ■  200  Gauss  (Inner  bullet),  (11)  Be 
•  50  fieuss  (dashed  outer  line). 
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Figure  4.  Effect  of  1.59  vertical  field  mismatch  In 
(a)  SLB,  (b)  Modified  betatron. 


Mismatch  In  Vertical  Field 


For  a  given  value  of  the  vertical  field  Bj,  the 
matched  Major  radius  of  the  bean  Is  r  •  svc/n2, 
where  r  Measures  the  distance  fron  the  Major  axis  of 
the  torus.  Transverse  oscillations  about  this  radius 
are  determined  fron  the  equation. 


4r  ♦  ar  (^  -  n,  iy  ♦  p2/d)  o|/y2 


where  b  is  the  Minor  radtus  of  the  drift -tube,  and 
ABz1s  the  arisaatch  In  the  vertical  field.  This 
equation  is  valid  in  the  thin-lens  Halt,  in  which 
the  solenoldal  lenses  act  as  a  continuous  radial 
focusing  force  on  the  beaa.  Equation  (2)  predicts  a 
Momentum  compaction  factor  •  ■  (1/2  -  n,  a*/b*  ♦ 
?/4)“*.  The  effect  of  a  1.59  vertical  field  Mis¬ 
match  Is  shown  in  the  particle  simulation  in  Fig.  4, 
In  which  we  assune  the  sample  UNK  beam  parameters 


given  above,  and  Ba  ■  100  Gauss.  For  these  para¬ 
meters,  «  ■  1,  so  that  the  bean  oscillates  about  an 
equilibrium  najor  radius  of  98.5  cm.  For  comparison. 


the  momentum  compaction  factor  of  the  weak-focuslng 
modified  betatron  Is  a  ■  (1/2  -  n«  a*/b*)*1  -  2.5. 
In  this  case  a  1.59  mismatch  in  the  vertical  field 


gives  the  beam  an  average  najor  radius  of  96  cm.  As 
seen  in  Fig.  4,  the  oscillation  about  this  position 
brings  the  bean  in  contact  with  the  wall.  Even  with 
the  Inproved  arisaatch  behavior  of  the  SLB,  however, 
keeping  the  vertical  field  Matched  to  the  beam  energy 
to  within  a  few  percent  will  require  delicate  tuning 
of  the  experiment. 


Stability  Properties 


A  circulating  high-current  electron  ring  nay  be 
subject  to  several  types  of  instabilities.  Including 
negative-mass,  resistive  wall,  and,  in  the  case  of 
the  SLB,  which  has  accelerating  gaps,  the  bean  break¬ 
up  instability.  In  addition,  single  particle  orbital 
resonances  nay  affect  bean  quality.  Here,  we  concen¬ 
trate  on  the  negative-mass  Instability,  since  It  Is 
potentially  the  most  serious  collective  instability.3 
An  analytic,  high-current  theory  of  this  and  some 
closely  related  Instabilities  is  described  In  a  com¬ 
panion  paper.*  The  theory  Is  directly  applicable 
only  to  devices  In  which  the  toroidal  coordinate  Is 
Ignorable  in  the  equilibrium,  such  as  the  modified 
betatron.  However,  we  use  the  theory  here  as  a  guide 
in  discussing  the  stability  of  the  SLB.  Our  numeri¬ 
cal  results  obtained  fron  3-D  simulations  using  IVORV 
are  not  restricted  this  respect.  However,  the  number 
of  simulations  we  can  perform,  and  their  length.  Is 
small  due  to  computing  expenses.  He  have  therefore 
concentrated  on  making  runs  at  high  currents,  where 
relatively  large  growth  rates  are  expected. 


The  first  case  we  look  at  is  a  10  kA,  5.5  MeV 
bean  with  a  •  2  cm,  and  R  •  1  n.  He  find  that  a  sol¬ 
enoldal  field  of  600  Gauss  can  transport  this  beam 
with  Minimal  (<109)  envelope  modulation.  The  solen¬ 
oldal  lenses  have  the  sane  configuration  as  In  Fig. 
2.  The  beam  is  given  a  small  initial  perturbation, 
and  the  growth  of  t  ■  20  fields  on  the  beam  Is  fol¬ 
lowed.  The  justifications  for  following  just  one 
node  In  the  simulations  are:  a)  In  the  linear  regime, 
node-coupling  should  be  negligible;  b)  By  following 
one  node,  direct  comparison  with  linear  theory  pre¬ 
dictions  for  that  particular  mode  are  possible;  c) 
Since  there  are  20  solenoids  around  the  torus,  the  t 
■  20  node  Is  likely  to  have  a  large  initial  perturba¬ 
tion;  d)  Computing  costs  severely  limit  the  number  of 
toroidal  nodes  one  can  afford  to  keep.  In  the  simu¬ 
lation  just  mentioned,  we  find  a  rapid  linear  growth 
rate  of  r  »  1.62* 10*  see-1.  The  effect  of  the  Insta¬ 
bility  on  the  beam  as  It  reaches  nonlinear  levels  Is 
shown  In  Fig.  5.  The  behavior  Is  similar  to  the  neg¬ 
ative-mass  behavior  that  one  would  expect  to  see  In  a 
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conventional  betatron  (Bg  >0).*  The  radial 
dtflactlon  of  the  bean  toon  In  Fig.  5(a)  It 
accompanied  by  tho  toroidal  bunching  In  Fig.  5(b). 
TMt  behavior  load*  ut  to  compare  tho  linear  growth 
rate  with  that  obtained  from  theory*  for  a  lg  •  0 
bean.  He  find  that  the  theoretical  growth  rate  It 
very  clote:  1.63k10*  tec-1.  In  Fig.  5.  we  compare 
result*  for  the  SLB  at  other  energies  with  the 
theory.  For  comparison,  results  for  the  tamo  bean  In 
a  Modified  betatron  with  i  1  tt  toroidal  field  are 
shown.  In  Fig.  7.  the  tame  curves  are  plotted  for  a 
1  kA  bean  with  a  •  1  cm.  For  these  parameters,  wa 
have  just  one  simulation  result,  at  5.5  MeV,  and  the 
growth  rate  Is  again  seen  to  be  close  to  the  ^  ■  0 
result. 
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are  plotted  In  (a). 
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Figure  6.  Growth  rates  of  I  •  20  mode  versus 
energy.  Simulation  results  for  SLB,  and 
analytic  results  for  conventional  (Bg  • 
0)  and  modified  (Bg  ■  1  kG)  betatronsare 
shown. 
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Figure  7.  Growth  rates  of  s  >  20  Instability  for  a  1 
kA  beam  with  1  cm  radius.  Curves  have 
same  meaning  as  In  Fig.  6. 


He  must  note  that  the  results  In  Fig.  €  for  the 
case  Bg  •  0  are  somewhat  artificial  when  y  <  18. 
In  this  regime,  weak  focusing  alone  Is  Insufficient 
to  hold  the  bean  together.  The  analytic  model  from 
which  the  growth  rates  are  obtained  uses  a  rigid  disk 
model  of  the  beam,  and  so  Ignores  the  force  balance 
required  within  the  beam.  It  Is  tempting  to  conclude 
from  the  results  In  Fig.  6  that  In  the  SLB,  the 
strong  focusing  provides  the  necessary  forces  the 
hold  the  beam  equilibrium  together,  but  that  the  neg¬ 
ative-mass  instabilities  on  the  beam  behave  as  If  the 
toroidal  field  were  taro.  Further  numerical  and  ana¬ 
lytic  work  will  be  performed  to  check  this  conclus¬ 
ion. 
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,\lle  have  developed,  a  new  linearized,  rigid-disk 

■Oder  of  negative  mass  Initabllltles  In  high  current 
betatrons.  Taking  the  beam  and  accelerator  cavity 
cross  sections  t.o  be  rectangular  permits  the  electro* 
magnetic  fields  to  be  evaluated  exactly  In  toroidal 
geometry.  Growth  rates  from  the  model  agree  well 
with  results  of  three-dimensional  numerical  simula¬ 
tions  for  beams  and  cavities  with  rectangular  or 
cylindrical  cross  sections.  Generally,  negative  mass 
Instability  growth  rates  are  greatest  for  beam  ener¬ 
gies  within  a  factor  of  two  of  the  so-called  transi¬ 
tion  energy  and  In  that  energy  regime  scale  Inversely 
with  the  square  root  of  the  toroidal  magnetic  field 
strength^  At  much  higher  energies,  growth  rates  are 
nearly  independent  of  toroidal  field  strength. 
Growth  rates  Increase  with  beam  current  and  toroidal 
mode  number,  but  the  scaling  laws  tend  to  be  compli¬ 
cated  by  competition  between  capacitive  and  Inductive 
components  of  the  toroidal  cavity  fields. 

'  Introduction  and  Summary 


Analytical  Results 

Me  consider  a  rectangular  cross  section  beam  and 
cavity  with  dimensions  as  shorn  In  Fig.  1.  The  beam 
may  be  located  radially  anywhere  In  the  cavity  but 
must  be  centered  axially  for  a  static  equilibrium. 
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Betatrons  employing  toroidal  magnetic  fields 
[1,2].  possibly  augmented  by  strong-focusing  fields 
[3-6],  seem  capable  of  accelerating  moderate  current 
electron  beams  to  high  energies.  Beam  stability, 
along  with  Injection  and  extraction.  Is  an  Important 
Issue  not  yet  completely  resolved.  Among  Instabili¬ 
ties,  the  negative  mass  mode  probably  Is  most  seri¬ 
ous,  as  It  has  been  for  other  electron  ring  devices 
[7].  In  this  paper  we  present  an  Improved  linear 
model  of  the  negative  mass  Instability  In  cold  beams, 
obtain  growth  rate  expressions,  and  compare  these 
findings  with  the  results  of  three-dimensional  com¬ 
puter  simulations. 

As  In  earlier  treatments,  we  approximate  the 
beam  as  a  string  of  rigid  disks  [8,9],  because  the 
Instability  has  the  form  of  a  uniform  transverse  beam 
displacement.  The  electromagnetic  fields,  however, 
are  evaluated  exactly,  although  for  a  beam  and  cavity 
of  rectangular  minor  cross  section  In  order  to  sim¬ 
plify  the  analysis.  Accurately  determining  the 
fields  Is  critical  because  of  delicate  cancellations 
among  terms  In  the  resulting  dispersion  relation 
[9].  Simulations  (with  IVORY  [10])  suggest  that  the 
choice  of  cross  section,  irfille  artificial,  does  not 
appreciably  alter  the  character  of  our  results. 

The  most  Important  outcome  of  the  new  model  Is 
Its  much  Improved  agreement  with  simulation  data, 
which  enhances  our  confidence  In  predictions  of 
experimental  behavior. 

Interestingly,  the  model  also  shows  that  Induct¬ 
ive  electric  fields  can  dominate  electrostatic  fields 
at  large  toroidal  mode  numbers,  giving  rise  to 
reduced  Instability  growth  rates  and,  for  a  few 
modes,  a  stability  window.  The  window  can  be  shifted 
to  low  toroidal  modes  by  moving  the  electron  beam 
toward  the  Inner  wall  of  the  cavity.  The  analogous 
situation  for  a  fast  rotating  e-layer  Immersed  In  a 
strong  azimuthal  magnetic  field  Is  treated  In  Ref. 
[11]  to  check  the  betatron  results  and  to  provide 
more  Insight  Into  the  underlying  physics. 


Figure  1.  Beam  and  toroidal  cavity  cross  sections. 

Treating  the  beam  as  a  set  of  rigid  disks  and 
evaluating  Its  fields  by  a  Green's  function  approach 
leads  to  a  dispersion  relation  of  the  usual  form 

B  2 

(q2  -  wz2)(a2  -  wr2  -  -Ji — •)  -  o2  -K  -  o  (i) 

0  -  C  T 

Where  Q  ■  w-tw0,  *■  •  V#/R,  R  •  (w1-fw2)/2,  and  Bg 
Is  the  polodial  cyclotron  frequency.  Expressions  for 
the  axial  and  radial  bounce  frequencies,  the  longi¬ 
tudinal  dielectric  function,  and  the  coupling  func¬ 
tion  are  extremely  complex  but  may  be  approximated 
as,  respectively,  [10] 

“z  "  "  *o2  ’  «1  ^7  <2> 

wp2  •  (1  -  n)  m2  -  92  -tV  (3) 

y  o 

2 

e  "  (93  J7  -  »4  “?) 

X  -  (yw0a  ♦  -^7  J7  %)2  "  tV  (°2  *  «)  (5) 

Here,  n  Is  the  betatron  field  Index,  v  Is  the  beam 

current  normalized  to  17  kA,  y  Is  Its  relativistic 
energy,  and  gj  are  geometrical  factors  of  order 
unity.  Much  of  the  physics  Is,  of  course,  contained 
In  those  factors,  the  latter  three  of  which  must  be 
accurately  evaluated  numerically  due  to  cancellations 
between  the  Inductive  and  electrostatic  components  of 
Eq.  (4)  and  (5). 
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Approximate  analytical  solutions  to  Eqs.  (I)-(5) 
can  bo  obtalnod  at  energies  above  the  transition 
energy,  1.  «..  for  the  right  side  of  Eq.  (2)  posi¬ 
tive.  See  Table  1.  It  Is  often  true  that 

*  *  J?  *0  <6) 

where  c*  Is  c  evaluated  at  Q  •  0.  When  this  Is  so. 
Cases  1A  and  1C  of  Table  1  correspond  to  the  usual 
negative  mass  growth  rate  expressions  for  a  low  cur¬ 
rent  betatron  without  [12]  and  with  [13]  a  toroidal 
magnetic  field.  The  pololdal  bounce  frequency  used 
In  the  table  Is  defined  as 

*B2  “  -r2  -z2/(Be2/^)  <7> 
TABLE  1.  APPROXIMATE  DISPERSION  RELATION  SOLUTIONS 
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Numerical  Results 


Me  next  solve  the  full  dispersion  relation, 
using  the  methods  described  In  Ref.  [10],  to  Illus¬ 
trate  typical  instability  growth  rate  scaling  with 
beam  current,  energy,  toroidal  field,  and  mode 
number.  We  choose  a  square  beam,  b  ■  1.76  cm,  cen¬ 
tered  In  a  square  cavity,  a  •  8.8  cm.  The  major 
radius  Is  R  ■  100  cm,  and  the  field  Index  Is  n  • 
1/2. 


Figures  2  and  3  depict  respectively  t  ■  1  and  20 
negative  mass  Instability  growth  rates  for  10  kA  and 
1  kA  beams  In  a  1  kG  guide  field.  Results  for  t  ■  1 
with  v  >  10  agree  well  with  Case  IB  of  Table  1.  The 
growth  rate  scales  linearly  with  t  and  /v.  Note  that 
the  growth  rates  at  high  energies  vary  roughly  as 
v’1/*  for  »  •  1  because  gj  slightly  exceeds  g*,  and 
as  y“5'2  for  t  •  10  (not  shown)  because  gj  and  g* 
are  for  that  case  equal.  For  high  energies  and  mod¬ 
erate  toroidal  magnetic  fields. 


Scaling  of  growth  rates  at  low  energies  Is  not  so 
clean  but  appears  to  go  as  (t/v)1/2. 

When  t  >  11  the  standard  negative  mass  Instabil¬ 
ity  Is  cut  off  by  Inductive  effects  for  x  negative, 
1.  e.  y  greater  than 

Teo  *  (*  -  *3'«4 )"1/Z  <»> 


The  cutoff  shown  for  the  1  kA  beam  at  t  ■  20  In  Fig. 
3  satisfies  Eq.  (9)  well.  Note  also  the  first 
appearance  of  a  second  Instability  branch  In  a  narrow 
band  around  y  ■  24  due  to  resonant  coupling  of  the 
beam  longitudinal  and  transverse  modes.  This  hybrid 
mode  Is  much  stronger  In  the  10  kA  beam  and  connects 
directly  to  the  usual  branch  of  the  instability  at  y 
■  18.  For  t  •  30  (not  shown),  yco  "  12;  the  hybrid 
negative  mass  Instability  Is  well  developed  for  both 
10  kA  and  1  kA,  and  Is  accurately  represented  by  Case 
2B  of  the  table. 


Figure  2.  Negative  mass  t  ■  1  growth  rates  for  10 
and  1  kA  beams  In  a  1  kG  guide  field. 


Figure  3.  Negative  mass  t  »  20  growth  rates  for  10 
and  1  kA  beams  In  a  1  kG  guide  field. 


Instability  growth  rate  scaling  with  Bg  Is 
presented  In  Fig.  4  and  5  for  the  10  kA  beam  at  t  *  1 
and  20.  Results  for  1  kG,  10  kG,  and  100  kG  guide 
fields  are  shown.  Corresponding  data  for  0  kG  are 
not  given  for  t  ■  1,  and  for  i  •  20  below  yco* 
because  they  are  Indistinguishable  from  the  1  kG 
results  at  energies  for  which  an  equilibrium  exists 
(above  about  18  MeV).  The  0  kG  betatron  Is  stable 
above  yco  for  A  *  20.  In  general,  we  expect  the 
guide  field  to  have  negligible  Influence  on  beam  sta¬ 
bility  for  low  t  (e.  g..  Fig.  4)  so  long  as  the 
pololdal  bounce  frequency  Is  less  than  ✓*.  For  the 
present  parameters  this  is  true  for  B»  lass  than  a 
few  kG  at  y  ■  10  and  a  few  thousand  kG  at  y  ■  50. 
Hence,  we  see  a  steady  drop  In  the  peak  of  the  growth 
rate  curve  as  Bg  Is  Increased  above  1  kG  until  It 
reaches  Case  1C  at  100  k6.  Note  that  growth  below 
the  transition  energy  Is  eliminated.  In  contrast, 
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••lew  rco  Instability  behavior  at  t  ■  20,  dis¬ 
played  In  Fig.  5,  is  qualitatively  similar  to  that  at 
•  »  1.  Growth  rates  above  the  Inductive  transition 
energy  are  expected,  based  on  Table  1,  to  fall  off 
with  Magnetic  field  as  V1,a  once  «g  somewhat 
exceeds  «a.  This  happens  above  a  few  kG  for  energies 
of  Interest  In  Fig.  5. 


Figure  4.  Negative  mass  t  •  1  growth  rates  for  a  10 
kA  bean  In  1,  10,  and  100  kG  guide  fields. 


Figure  5.  Negative  mss  t  •  20  growth  rates  for  a  10 
kA  baa*  In  1,  10,  and  100  kG  guide  fields. 


Figure  4.  Negative  mss  s  •  1  growth  rates  for  a  10 
kA  bean  in  a  1  kG  guide  field. 


Comparison  with  Slaxilatlons 

No  conclude  this  paper  with  a  brlaf  comparison 
of  dispersion  relation  curves  with  growth  rates 
determined  by  the  computer  simulation  code  IVORY. 
The  simulation  data  points  were,  for  the  most  part, 
already  reported  in  Ref.  [9].  Figure  6  gives  s  ■  1 
results  for  a  10  kA  bean  In  a  1  kG  guide  field. 
Agreement  Is  excellent  except  for  the  r  *  5  data 
point.  The  discrepancy  probably  Is  due  to  beam  temp¬ 
erature  In  the  simulation.  Figure  7  treats  the  same 
beam  but  with  y  held  fixed  at  12  and  s  varied  from  1 
to  20.  Again,  agreement  between  theory  and  simula¬ 
tion  Is  excellent.  The  few  simulations  in  Fig.  7 
performed  with  rigid  disks  Instead  of  discrete  par¬ 
ticles  further  vindicate  the  rigid  disk  approximation 
made  In  our  model. 


Figure  7.  Negative  mss  growth  rates  at  various  s 
for  a  10  kA  beam  in  a  1  kG  guide  field. 
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